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Abstract 

The NCC code was validated for a case involving staged transverse 
injection into Mach 2 flow behind a rearward facing step. Comparisons 
with experimental data and with solutions from the FPVortex code were 
performed to assess the performance of the NCC code. The code was then 
used to perform computations to study fuel-air mixing for the combustor 
of a candidate rocket-based combined cycle engine geometry. Comparisons 
with a one-dimensional analysis and a three-dimensional code (VULCAN) 
were performed to assess the qualitative and quantitative performance of 
the NCC solver. 

1. Introduction 

NASA has been focusing on several alternative propulsion concepts in 
order to reduce the cost of access to space. Under NASA’s Integrated Space 
Plan three key technologies are being studied: scramjet, Rocket-Based 

Combined Cycle (RBCC) and Turbine-Based Combined Cycle (TBCC). 

The Integrated Systems Test of an Airbreathing Rocket (ISTAR) 
project is proposing to use the RBCC concept of combining rocket-based 
propulsion with dual-mode RAM/SCRAM air-breathing propulsion to cre- 
ate a strut-jet engine. This engine would combine rocket-containing struts 
with conventional combustor ducts to obtain high combustion efficiencies 
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in all flight regimes. RBCC engines take advantage of the synergistic in- 
teractions between the rocket and airbreathing elements of the engine and 
the use of high specific-impulse cycles to yield a mission-averaged specific 
impulse that is higher than all-rocket technology. 

High-fidelity predictions of the physics for RBCC-combustor phenom- 
ena could be attained by the solution of Navier-Stokes equations for reacting 
flows, and this is a very challenging computational task. The complexity of 
fuel-air mixing in conjunction with large heat-release due to combustion is 
magnified by the combination of high-speed rocket exhaust flows and the 
relatively low-speed duct flows. 

Recent work by Steffen et. al 1 has focused on the analysis of a typical 
RBCC system during rocket-only operation. The authors had performed an 
extensive numerical study to determine the parameters that would affect 
rocket-mode RBCC performance. The authors had determined that the 
expansion and combustion process downstream of the rocket nozzle could 
greatly impact the overall cycle specific-impulse. Stewart and Suresh 2 have 
performed multidisciplinary analysis of an ISTAR/RBCC configuration us- 
ing the VULCAN code 3 . Their results will serve as a basis of comparison 
for the results in the current report. Daines and Segal 4 have identified 
rocket-enhanced mode performance, flame-holding and thermal choking as 
important issues for RBCC engines. These issues can only be effectively 
studied by detailed, high-fidelity analyses of RBCC combustor dynamics. 

One of the challenges facing the designers in the ISTAR project is 
the lack of tools to quickly assess the combustion efficiency of proposed 
designs. In a typical design cycle, the design team may typically require 
answers to questions like - what impact will a particular fuel/air ratio have 
on combustion efficiency and how will the location and size of the injectors 
affect mixing in the combustor? 

These and many other questions require robust and accurate compu- 
tational tools which can provide reliable answers for the designers in order 
to help them reduce the size of the design space and also eliminate ‘show- 
stopping’ design paths early in the design cycle. The purpose of this report 
is to document the performance of a combustion simulation tool, the Na- 
tional Combustion Code (NCC), in providing answers to combustion design 
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issues for RBCC-type applications. 

The report is organized into four sections. The introduction section is 
followed by a very brief description of the NCC code — the governing equa- 
tions, spatial discretizations, and turbulence and chemistry modeling. This 
section is followed by validation studies for staged transverse injection into 
a supersonic throughflow behind a rearward-facing step. A comparison of 
computational and experimental results are then presented for the valida- 
tion case. The computational results for a candidate RBCC-combustor are 
then presented. The results are compared with the computational results 
of two other codes. The final section is a summary of the present work. 

2. Computational Tools 

The National Combustion Code (NCC) 5 is a state-of-the-art compu- 
tational tool which is capable of solving the time-dependent, Navier-Stokes 
equations with chemical reactions. The NCC is being developed primarily 
at NASA Glenn in order to support combustion simulations for a wide range 
of applications, and has been extensively validated and tested for low-speed 
chemically reacting flows. 


2.1 Governing Equations 


The Navier-Stokes equations for unsteady, compressible, chemically re- 
acting gas flow with N species can be written as 


dQ dQ d(F-F v ) 

dr dt dx 

d(G-G v ) d(H-H v ) = 

+ dy + dz 

where the vectors Q, F,G, H, F V ,G V , H v and / are defined as 

Q = (p, pu, p% pw , pE, pK 1 pe, pY \ , ■ ■ ■ , pY N _ i) T 

F = (pu, pu 2 + p g . puv , puw, ( pE + p)u, 

\ T 1 

pun ? pue ? puYi , ■ ■ ■ ? puYn- i) 
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I=(0,S u ,S v ,S m ,S h ,S li ,S e ,S 1 ,%S i ,;S N . 1 ) T 

where p,p g ,u,v,w, E, K,e and E*- represent the density, gauge pressure, x-, 
y-, z- cartesian velocity components, total energy (e + \(u 2 + r> 2 + w 2 )), 
turbulent kinetic energy (TKE), dissipation rate of TKE, and species mass- 
fraction, respectively. The source terms S u , S v , S w , Sh, S K , S e arise from a 
non-linear k — e turbulence model, and the source terms Si are the chemical 
source terms for species i. 


In order to overcome potential convergence problems for low Mach 

number flows, a preconditioning term (T^y ) has been added to the standard 
Navier-Stokes equations 6 . During the solution procedure, this term is driven 
to zero as the solver marches in pseudo- time, r. In addition, the pressure 
term has been decomposed as a reference pressure (pQ ) and a gauge (or 
fluctuating) pressure (p g ) . i.e. p = po -\- p g ■ The vector Q is thus the vector 
of primitive dependent variables 


Q = (p 


9- 


U. V , w, 


K 


«,e,y l 
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and the preconditioning matrix T is written as 


- 1 /(3 0 0 

u/j3 p 0 

v/ (3 0 p 

w/(3 0 0 

H/ j3 — 1 pu pv 

k//3 0 0 

e/(3 0 0 

Y\f (3 0 0 

Y 2 /(3 0 0 


lY N -il(3 0 0 
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0 . . p . . . 0 0 
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where (3 is a scaling factor for the eigenvalues, and is computed as (3 = 
u 2 + v 2 + w 2 . This ensures that all eigenvalues are of the same order of 
magnitude 7 . 


2.2 Time-Integration 


Equation (1) is integrated volumetrically by applying a cell-centered, 
finite- volume discretization. 


r $Q A , d Q A , 

r — ^ Av + — -An+ 

OT Ot 

nf nf 

j2(F-F v )A x + J2(G 

n — 1 n — 1 

nf 

H = IAv 

n — 1 


G V )Ay + 


(2) 


where Av is the volume of the computational element, and nf is the number 
of faces constituting the element. A x , A y , A z are the projections of the 
face area in the x, y, and z directions, respectively. Second-order accurate 
central-differences are used for the inviscid and viscous flux discretizations, 
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and a Jameson operator 8 (a blend of 2nd and 4th-order dissipation terms) 
is used to maintain numerical stability. In order to enhance convergence 
acceleration in pseudo-time, implicit residual smoothing is used to smooth 
the computed residuals. Turbulence-closure is obtained by a low- Reynolds 
number two-equation k — e model 9 . For a steady state solution, ^ is 
set to zero, and equation (1) is integrated in pseudo-time with an explicit 
four-stage Runge-Kutta algorithm. 


3. Test Results and Discussion 

The National Combustion Code (NCC) is used to perform simulations 
of cross-flow injection in high-speed flows. A validation case was used to 
establish the performance of the NCC and the results compared with exper- 
imental and computational data. A candidate RBCC-combustor geometry 
was then computed with the NCC and the results compared with two other 
computational codes. 


3.1 Validations Study: Staged Transverse Injection 

The staged normal injection of two jets located behind a rearward- 
facing step into a Mach 2 freestream airflow is computed as the validation 
case. Experimental data for this configuration was obtained by McDaniel 10 . 
The computational setup for this problem is shown in figure 1. The mesh, 
with about 150,000 tetrahedal elements, is identical to the mesh used to 
compute this problem with the FPVortex code 11 . This would allow a one-to- 
one comparison between the current NCC computation and the FPVortex 
computation. 

The geometry consists of a channel with a 3.18mm(=H) step and two 
1.93mm(=D) injectors positioned on the channel centerline at 3 and 7 step- 
heights, respectively, downstream of the step. The channel height is 11.03 
injector diameters (D) and the width is 15.79 injector diameters. The 
freestream flow conditions correspond to a stagnation temperature of 300K 
and stagnation pressure of 274KPa. The inflow boundaries for the main 
flow and the two injectors were set as fixed, specified inflow conditions. All 
flow variables were extrapolated at the outflow boundary. 
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The NCC computed contours on the z=0 plane (centerline of jets) for 
axial velocity and static pressure are shown in figure 2. The penetration of 
the two jets into the main flow is indicated in figure 3 (marker concentration 
contours of the two jets). Mach number contours and the concentration of 
the main flow marker are shown in figure 4. 

Comparisons of the experimental data with the NCC and FPVortex 
computed normalized streamwise velocity and normalized static pressure 
profiles at various axial locations are shown in figures 5- 14. The loca- 
tions of the upstream and downstream jets correspond to X/D=0.0 and 
X/D=6.6, respectively. The normalized velocity and pressure predicted by 
the NCC is qualitatively similar to the predictions of the FPVortex code, 
with some differences in the near wall regions. Note that the computa- 
tional grid and flow and boundary conditions are identical for the NCC 
and FPVortex codes. The only probable difference between the two com- 
putations is the turbulence model. Overall, it may be concluded that the 
current NCC computations compare fairly well with available experimental 
and computational data. 


3.2 RBCC Combustor Computations 

Scram combustion within an RBCC-combustor configuration is mod- 
eled as a test case to evaluate the computational performance of the NCC 
code for high-speed combustor flows. The ‘E-201’ combustor configuration 
of the ISTAR vehicle is modeled, and the results are qualitatively compared 
with the computations performed with the VULCAN code 3 by Stewart and 
Suresh 2 . Comparisons are also made for certain flow quantities with a 1- 
D analysis with the FAST code 12 . The current NCC computation focuses 
on studying the mixing of the fuel (normal injection) with the secondary- 
air/rocket- exhaust main-stream flow. 

The schematic for the current NCC computations is shown in figure 15. 
The current computation corresponds to a flight Mach Number of 3.5 (50K 
ft) with p tot —77 Apsi and T tot — 1345. 5R. These conditions were used to spec- 
ify the inflow conditions for the secondary flow (air). Subsonic flow con- 
ditions with p tot and T tot held fixed and Mach number allowed to float 
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(with an initial value of M=0.63) were used. The flow conditions for the 
(six) rocket faces correspond to rocket exhaust conditions of u=7930ft/s, 
T=2390R, p=0.02391bm//f 3 . The rocket faces were specified with fixed, 
supersonic inflow conditions. The specified mass fractions of the species 
are C 2 H 4 = 0.05162, 0 2 = 0.15885, C0 2 = 0.2289, H 2 0 = 0.56059. The 
chemical reactions were modeled with a 6-species, 3-step, finite-rate reduced 
mechanism for ethylene described below. No turbulence-chemistry interac- 
tions were modeled in the current simulation. 


Reactants 

Af 

°7 

Ef 

C 2 H 4 + 0 2 ^ 2 CO + 2 H 2 

2.10E14 

0.00 

1.80153E4 

2 CO + 0 2 ^ 2 C0 2 

3.48E11 

2.00 

1.013486E4 

2 H 2 + 0 2 ^ 2 H 2 0 

3.00E20 

-1.00 

0.0 


The fuel is injected normal to the main flow (mixture of secondary 
air and rocket exhaust) at a downstream location through two transverse 
injectors of 0.331in effective diameter. The injection is sonic, with p tot = 
531.2psi, T tot = 1260 R and Ui n j = 1615.22ft/ s. The injectors are gridded 
as triangular slots with matching area, momentum, and mass flow. 

Computational results for the NCC code were compared with two pre- 
vious computations with a 1-D code (FAST) and a 3-D code (VULCAN) 
as reported by Stewart and Suresh[2]. The VULCAN simulation used 1.02 
million nodes, and was spatially first-order accurate. The current NCC 
computation uses 85,154 elements and is spatially second-order accurate. 
The chemistry model and flow conditions used for the current computation 
are identical to those used in the VULCAN simulation. A tabular compari- 
son of area-averaged secondary-inflow quantities for the three computations 
is shown in Table 1. All the quantities are non-dimensionalized with the 
respective values reported from computations with the FAST code. 
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Quantity 

FAST 

VULCAN 

NCC (current) 

mdot 

1.00 

1.170 

1.096 

ptot 

1.00 

1.0373(fixed) 

1.0373(fixed) 

Ttot 

1.00 

1.0 193 (fixed) 

1.0 193 (fixed) 

Mach Number 

1.00 

1.2365 

1.0635 

Static Pres 

1.00 

0.9073 

0.9642 

Static Temp 

1.00 

0.9789 

0.9837 


Table 1 . Area- averaged secondary-inflow quantities 
A tabular comparison of area-averaged flow quantities at the internal 


nozzle exit plane is 

shown in 

Table 2. 



Quantity 

FAST 

VULCAN 

NCC (current) 

mdot 

1.00 

1.1239 

1.0264 


Mach Number 

1.00 

0.9546 

0.9897 


Static Pres 

1.00 

1.2426 

1.6800 


Static Temp 

1.00 

0.9896 

1.0094 



Table 2. Area- averaged flow quantities at nozzle exit plane 

A plot of the static temperature contours in a plane containing the 
centerline of the injector face(s) as computed by Stewart and Suresh with 
the VULCAN code is shown in figure 16. The plot shows that there is 
considerable penetration of the fuel jets into the main chamber, and fairly 
uniform burning downstream of the injectors. Note that the fuel penetration 
reported in many experiments (such as the one shown in figure 3) is typically 
more limited than the penetration shown in figure 16. 

Figure 17 shows the fuel mass-fraction and temperature contours for 
the current NCC computation at y/Y=0.5 (injector center-line plane). Note 
that ‘Y’ is the distance between the symmetry face and the side-wall (see 
fig. 15). Also, y/Y=0 at the symmetry face and y/Y=l at the side-wall. 
The current results differ from the Vulcan solution in two ways - the limited 
penetration of the fuel jets into the main flow, and the non-uniform burning 
downstream of the injectors. The current results show a distinct, low- 
temperature core-flow surrounded by two streams of hot flow in the vicinity 
of the upper and lower walls. The maximum temperature is 2247K, as 


NAS A/TM— 2005-2 1 3647 


9 



compared to 2870K reported in the VULCAN simulation. 

One major difference between the VULCAN and NCC simulation is 
the mesh size and the mesh distribution. The VULCAN simulation used 
a structured, multiple-block approach with 1.02 million mesh points while 
the current NCC computation uses a unstructured mesh approach with only 
85154 mesh points. This could possibly account for the differences in the 
two solutions. In order to somewhat further study this possibility, a second 
NCC computation was performed with a 225,204 element tetrahedral mesh. 
This finer mesh was obtained by doubling the number of points which rep- 
resent the edges of the computational domain. The finer mesh does provide 
somewhat improved near- wall resolution. The fuel mass-fraction and tem- 
perature contours for the finer grid computations are shown in figure 18. 
There is no significant difference between the coarse and fine grid solutions. 
The jet penetration and fuel-air mixing characteristics are qualitatively the 
same for the fine and coarse grid computations. Hence, the grid size and/or 
distribution is ruled out as the major source of the difference between the 
current NCC computations and the previously reported VULCAN results. 

Figures 19 and 20 show contours for fuel and H 2 O mass-fractions at 
four different spanwise locations of the RBCC combustor. These locations 
were chosen to study the spanwise variation (y-axis of figure 15) in the com- 
bustor flow- field. Recall that y/Y=0 at the symmetry face and y/Y=l at 
the side wall (fig. 15). The figures show that there is considerable variation 
in fuel-penetration and product formation at the different spanwise loca- 
tions. The lack of mixing between the flow near the walls and the center of 
the combustor is very apparent at all four locations. This would suggest that 
a better mixing strategy would probably be required to enhance mixing and 
combustion efficiency. Figures 21 and 22 show contours for CO and CO 2 
mass-fractions. The production of CO and CO 2 also varies considerably 
across the combustor. Figures 23 and 24 show contours for mach- number 
and pressure. The mach- number and pressure both change considerably 
as the flow moves downstream from the fuel injectors. The computed re- 
sults are further illustrated in figures 25-32, which show the distribution 
of static temperature, static pressure, axial velocity, Mach number, and 
various species fractions at different axial locations. 
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4. Conclusions and Future Work 


A successful computational study was conducted to validate and then 
apply the National Combustion Code (NCC) to high speed combustor prob- 
lems with transverse fuel injection. The RBCC-combustor computations 
showed that the NCC code could prove to be a viable tool for high-fidelity 
analysis of combustor geometries. The comparison of NCC computations 
with experimental data and three other computation codes demonstrated 
that the NCC code compared fairly well for the class of problems studied. 
While there were some significant differences in the detailed flowfield be- 
tween the NCC and Vulcan results for the RBCC-combustor geometry, the 
comparison of area-averaged quantities was fairly close. It would be inter- 
esting to revisit the comparisons when experimental data for this particular 
geometry becomes available. 
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Figure 1. — Computational setup for transverse staged (sonic) injection. 




u -86 -38 10 58 106 154 202 250 298 346 394 442 490 538 586 


(a) 


Axial velocity contours, m/s 



i ii i i 

p 13 489 33 704 53 918 74 132 94 347 114 561 134 776 154 990 

23 596 43 811 64 025 84 240 104 454 124 668 144 883 

Static pressure, Pa 
(b) 

Figure 2. — NCC computed contours on plane containing the jet centerline, for staged transverse injection, (a) Axial 
velocity, (b) Static pressure. 
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(a) 
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Marker concentration, upstream jet 


0.81 0.88 0.94 




0.06 0.12 0.19 0.25 0.31 0.38 0.44 0.50 0.56 0.62 0.69 0.75 0.81 0.88 0.94 


Marker concentration, downstream jet 
(b) 

Figure 3. — NCC computed contours of marker concentrations indicating penetration distance into the 
main flow, for staged transverse injection, (a) Upstream jet. (b) Downstream jet. 




0.06 0.12 0.19 0.25 0.31 0.38 0.44 0.50 0.56 0.62 0.69 0.75 0.81 0.88 0.94 


Marker concentration, main flow 
(b) 

Figure 4. — NCC computed contours of the main flow, for staged transverse injection, (a) Mach number, 
(b) Marker concentration. 
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Figure 5. — Comparison of upstream injector normalized Figure 6. — Comparison of upstream injector normalized 

velocity data, for staged transverse injection. x/D = 0.0. pressure data, for staged transverse injection. 

x/D = 0.0. 




Figure 7. — Comparison of upstream injector normalized Figure 8. — Comparison of upstream injector normalized 

velocity data, for staged transverse injection. x/D = 3.0. pressure data, for staged transverse injection. 

x/D = 3.0. 
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Figure 9. — Comparison of upstream injector normalized 
velocity data, for staged transverse injection. x/D = 6.0. 



Figure 10. — Comparison of upstream injector normalized 
pressure data, for staged transverse injection. 
x/D = 6.0. 




Figure 11. — Comparison of downstream injector normal- 
ized velocity data, for staged transverse injection. 
x/D = 6.6. 


Figure 12— Comparison of downstream injector normal- 
ized pressure data, for staged transverse injection. 
x/D = 6.6. 
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Figure 13. — Comparison of downstream injector normal- 
ized velocity data, for staged transverse injection. 
x/D = 12.8. 
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Figure 14. — Comparison of downstream injector normal- 
ized pressure data, for staged transverse injection. 
x/D = 12.8. 



Figure 15. — Computational grid and schematic for RBCC combustor. 
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Figure 16. — VULCAN code temperature contour in a plane containing the injector 
centerline (Ref. [2]). 
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Figure 17. — NCC code coarse grid solution for RBCC combustor in a plane containing the injector centerline. 


NAS A/TM— 2005-2 1 3647 


18 






Fuel 

>s fraction 

0.967 

0.901 

0.835 

0.768 

0.702 

0.636 

0.570 

0.503 

0.437 

0.371 

0.305 

0.239 

0.172 

0.106 

0.040 


Temperature, 
K 

2297.70 
2172.79 
2047.88 
1922.97 
1798.06 
1673.15 
1548.23 
1423.32 
1298.41 
1173.50 
1048.59 
923.68 
798.77 
673.86 
548.95 

Figure 18. — NCC code fine grid solution for RBCC combustor in a plane containing the injector centerline. 
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Figure 19. — NCC code fuel mass fraction for RBCC combustor at 4 spanwise locations. 

(a) Centerline of rocket plane (y/Y = 0.06). (b) Centerline of fuel injection plane (y/Y = 0.60). 
(c) Between rocket plane and fuel injection plane (y/Y = 0.36). (d) Between fuel injection 
plane and side wall (y/Y = 0.79). 
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Figure 20. — NCC code H 2 0 mass fraction for RBCC combustor at 4 spanwise locations. 

(a) Centerline of rocket plane (y/Y = 0.06). (b) Centerline of fuel injection plane (y/Y = 0.60). 
(c) Between rocket plane and fuel injection plane (y/Y = 0.36). (d) Between fuel injection 
plane and side wall (y/Y = 0.79). 
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Figure 21 . — NCC code CO mass fraction for RBCC combustor at 4 spanwise locations. 

(a) Centerline of rocket plane (y/Y = 0.06). (b) Centerline of fuel injection plane (y/Y = 0.60). 
(c) Between rocket plane and fuel injection plane (y/Y = 0.36). (d) Between fuel injection 
plane and side wall (y/Y = 0.79). 
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Figure 22. — NCC code C0 2 mass fraction for RBCC combustor at 4 spanwise locations. 

(a) Centerline of rocket plane (y/Y = 0.06). (b) Centerline of fuel injection plane (y/Y = 0.60). 
(c) Between rocket plane and fuel injection plane (y/Y = 0.36). (d) Between fuel injection 
plane and side wall (y/Y = 0.79). 
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Figure 23. — NCC code pressure contours for RBCC combustor at 4 spanwise locations. 

(a) Centerline of rocket plane (y/Y = 0.06). (b) Centerline of fuel injection plane (y/Y = 0.60). 
(c) Between rocket plane and fuel injection plane (y/Y = 0.36). (d) Between fuel injection 
plane and side wall (y/Y = 0.79). 
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Figure 24. — NCC code Mach number contours for RBCC combustor at 4 spanwise locations, 
(a) Centerline of rocket plane (y/Y = 0.06). (b) Centerline of fuel injection plane (y/Y = 0.60). 
(c) Between rocket plane and fuel injection plane (y/Y = 0.36). (d) Between fuel injection 
plane and side wall (y/Y = 0.79). 
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Figure 25. — NCC computed temperature contours for RBCC combustor at axial locations. 



Figure 26. — NCC computed pressure contours for RBCC combustor at axial locations. 
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Figure 27. — NCC computed axial velocity contours for RBCC combustor at axial locations. 



Figure 28. — NCC computed Mach number contours for RBCC combustor at axial locations. 


NAS A/TM— 2005-2 1 3647 


27 



Figure 29. — NCC computed fuel mass fraction contours for RBCC combustor at axial locations. 



Figure 30. — NCC computed O 2 mass fraction contours for RBCC combustor at axial locations. 
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Figure 31 . — NCC computed C0 2 mass fraction contours for RBCC combustor at axial locations. 
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Figure 32. — NCC computed H 2 0 mass fraction contours for RBCC combustor at axial locations. 
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